Influence of absorbing dielectric background on bistable response of dense collection 

of two-level atoms 
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The stationary problem ol light interaction with a dense collection of two-level atoms embedded 
in a dielectric is treated semiclassically. The effect of dielectric absorption on the phenomenon 
of intrinsic optical bistability is discussed in terms of the complex local-field enhancement factor. 
The calculations accounting propagation effects (nonlocal regime) show that even low dielectric 
absorption results in elimination of response instabilities. 
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I. INTRODUCTION 

The effect of optical bistability is the corner-stone of 
optical information processing. It is the basis for imple- 
mentation of such devices as all-optical switches, tran- 
sistors, and logical gates. The straightforward method 
to obtain the bistable response is to place a nonlinear 
medium inside a feedback system, for example, Fabry- 
Perot resonator [l|. In the late 1970s Bowden and Sung 
introduced [|| the idea of the mirrorless (intrinsic) opti- 
cal bistability as a result of interatomic (dipole-dipole) 
correlations. This idea was then developed in a number 
of papers (see, for example, 

The usual starting point in the analysis of intrin- 
sic optical bistability (IOB) is the semiclassical gener- 
alized Bloch equations [6| which take into account the 
near-dipole-dipole interactions of two-level atoms (the so- 
called local field correction). However, there is a problem 
of influence of the background (host) dielectric medium 
on the bistable response. Crenshaw proved [7[ that the 
semiclassical approach gives the correct explanation of 
local-field enhancement by the linear dielectric in ac- 
cordance with the microscopic quantum electrodynam- 
ics (QED). Interestingly, this enhancement cannot be ac- 
counted for by the macroscopic QED. Therefore, if the 
process of spontaneous emission is not to be examined, 
we can use the semiclassical Bloch equations. 

The local-field enhancement due to the host dielectric 
is measured by the factor £ = (e + 2)/3, where e is the 
dielectric permittivity. Since e can be complex, we can 
automatically take into account dielectric absorption. It 
was shown in Ref. Q that the real part of £ reduce 
the IOB threshold while the imaginary part acts in the 
opposite direction. In one of our previous publications 
[9| we studied the effects of local and nonlocal bistability 
and internal coherent reflection under the assumption of 
absence of the background dielectric influence. In this 
short paper we try to fill this omission. 

The paper is divided into several parts. Section II 



is devoted to the solution of the stationary problem of 
monochromatic radiation interaction with the dense col- 
lection of two-level atoms in a dielectric. We analyze the 
influence of dielectric properties on the main parameters 
of local bistability. In Section III we consider the overall 
permittivity of such a media and its bistable response. 
Finally, in Section IV the problem of radiation propaga- 
tion in a layer of finite thickness is discussed. The effect 
of dielectric on the so-called nonlocal bistability is nu- 
merically analyzed. 



II. BISTABLE RESPONSE 

In the paper we consider the normal incidence of lin- 
early polarized monochromatic plane wave on a dense res- 
onant medium. We start with semiclassical Bloch equa- 
tions which take into account near dipole-dipole (NDD) 
interactions of two-level atoms as well as background di- 
electric effect @,H, 
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where W is the inversion (population difference of two 
levels), R is the atomic polarization, E is the amplitude 
of macroscopic electric field, fi is the transition dipole 
moment, A is the detuning of radiation from resonance, 
71 and 72 are the population and polarization relaxation 
rates respectively, W eq is the value of inversion at equi- 
librium. The parameter 
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is responsible for NDD interaction strength, N is the den- 
sity of two-level atoms per unit volume, % is the Planck 
constant. Here £ — (e + 2)/3 is the local-field enhance- 
ment factor due to polarizability of host material with 
dielectric constant e (generally, complex). Appearance 
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FIG. 1: Dependence of inversion on the light amplitude at 
(a) real values of factor £, (b) complex values of factor £. 
Calculation parameters: W cq = — 1, S — —2, b — 4, T — 0.1. 
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FIG. 2: (Color online) Positions of the inflection points as a 
function of the NDD parameter b. Other parameters: W eq = 
-1, 5 = -2. 
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FIG. 3: (Color online) The region of bistability existence in 
the plane (S, b). 
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FIG. 4: The dependence of minimal NDD parameter value 
(a) and corresponding detuning (b) on the imaginary part of 
the local-field enhancement factor £. 



hancemcnt of the magnitude and the phase shift (if e is 
complex) of the electric field; (ii) Lorentz frequency shift 
(u>L = tW) enhancement due to real part of £; (iii) coop- 
erative decay due to imaginary part of £ in both equations 
and ©. 

In stationary regime we obtain the expression for the 
polarization, which can be written in dimensionless no- 
tation as follows 



R = 



1 



awn 



2i(5-£bW) - 1 



(4) 



where f2 = fiE/firfz, S = A/72, b = 6/72- It is seen that 
the polarization depends on the inversion which can be 
calculated from the equation 



of factor £ results in three local-field effects j8j: (i) en- 



T(W - Weq) \i(S - £bW) - 1| 2 + W\£Vl\ 2 = 0. (5) 
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Here V = 71/72- 

The cubic equation describes the phenomenon of 
intrinsic optical bistability (IOB) well known from the 
previous researches. However, the presence of dielec- 
tric background changes the quantitative characteristics 
and conditions of IOB. As can be seen from Fig. HJa), 
the Lorentz shift enhancement (the second effect named 
above) results in easier observation of bistable response. 
In vacuum case {£ = 1), if the other parameters are the 
same, there is no bistability at all. On the other hand, 
if the dielectric is absorptive [complex factor £, see Fig. 
[TJb)], bistable loop not only shifts towards higher light 
intensities, but also gets more and more narrow. 

This implies that absorption of the background mate- 
rial changes the region of IOB existence somehow. Let 
us consider this question in detail. The necessary con- 
dition for the Eq. §5§ to demonstrate bistability is that 
the equation <i|$l| 2 /G?VT = has two different roots in 
the physically appropriate range — 1 < W < 0. This 
equation can be written in explicit form as 



2b 2 \e\ 2 W 3 - b (2lm£ + 28Re£ -\ 
+ W cq (5 2 + 1) = 0. 



! W cq ) W 1 



(6) 



Two solutions of Eq. ^ in the range — 1 < W < rep- 
resent the inflection points of bistability which are seen 
in the Fig. Q](tlie points of jump from one solution to an- 
other). When these points coincide, the bistability loop 
disappears. Figure [5] shows the behavior of the inflec- 
tion points as NDD parameter is decreasing. The critical 
(minimal) value of b at which bistability still exists de- 
creases for larger real parts of the local-field enhancement 
factor and increases for larger imaginary parts of I. At 
the same time the inversion corresponding to the jumps 
between the solutions is growing. 

This means that the region of parameters S and 6, 
where the bistability, occurs changes as well. Some exam- 
ples of these regions of bistability existence are demon- 
strated in Fig. [3J above the curves plotted the bistabil- 
ity is present. In the case of vacuum {£ = 1) we obtain 
the curve that was reported previously Q. The minimal 
value of the NDD parameter in this case is fr m i„ = 4 for 
the detuning <$ min = — 1. If £ increases and is still real, 
the region of bistability is getting wider. For example, 
for t — 2 we have 6 m j n = 2 at the same value of the de- 
tuning (5 m in = — 1. However, if the background dielectric 
is absorptive (I is complex), this minimal value of NDD 
parameter increases. 

From Fig. [3] one can see that & m ; n for lm£ = 1 is larger 
than for Im£ = 0. However, for Imf = 2 it is already 
smaller than in the former case. Detailed study of the 
behavior of the minimal value of the NDD parameter 
[Fig. |U(a)] shows that there is a maximum in dependence 
6 m i n (Im^). As Im£ tends to infinity, the b mm is slowly 
decreasing down to zero. At the same time the value of 
detuning, corresponding to the minimal NDD parameter, 
monotonically tends to zero, too. Note that, for smaller 
Re£, the curve & m i n (Inu?) has more pronounced maximum 
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FIG. 5: Dependence of (a) real and (b) imaginary parts of 
the dielectric permittivity on the light amplitude. Calculation 
parameters: W cq = —1, 5 — —1, b — 4, V = 0.1, Rel = 1.5, 
Im£ is a varying parameter. 



and then slower converge to zero. In contrast to this, 
<5 m in — > faster for smaller Re£. 

Moreover, as lm£ is growing, the bistability existence 
region, at first, gets much more narrow for 5 < S mlu and 
slightly wider for S > <5 m j n (see Fig. [3]). But then it tends 
to get wider in both cases. In the limit lm£ —> 00, the 
region of bistability occupies all the half-plane b > 0. 



III. DIELECTRIC PERMITTIVITY 

The resulting dielectric permittivity of the dielectric 
doped with two-level atoms can be written as 



Etotal = 1 + 47TX, 



(7) 



where \ = P/E is the susceptibility, P is the total polar- 
ization defined as a sum of the linear polarization of the 
background dielectric (permittivity e) and the nonlinear 
polarization P rcs = 2NfiR due to resonant atoms, so that 
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FIG. 6: Dependence of reflection at the vacuum-medium in- 
terface on the light amplitude. Calculation parameters are 
the same as in Fig. [5] 
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FIG. 7: Dependence of reflection (a, b, c) and transmission (d) 
on the light amplitude. Calculation parameters: Wcq = —1, 
5 = -1, 6 = 4, T = 0.1, light wavelength A = 0.5 fim. The 
enhancement factor (a) I = 1, (b, d) £ = 2, (c) I — 3. The 
thickness of the layer is L — 0.5 /mi. 



(8) 



Substituting the stationary expression (0} to the equa- 
tions ([8]) and ([7]) we obtain 



£total 



3il 2 bW 



i{S-£bW) - 1' 



(9) 



It is worth to recall that I — (e + 2)/3. 

Expression ([9]) with cubic equation (0 defines the 
bistable behavior of the dielectric permittivity (Fig. [5]). 
This bistability, in accordance with the results of the pre- 
vious section, is getting more and more narrow at first 
and then wider, as the imaginary part of local-field en- 
hancement factor is growing. Moreover, it is seen that 
the real part of Etotai appears to be negative at low- 
intensive branch of the bistability. This negativity cor- 
responds to the well-studied effect of the internal coher- 
ent reflection [91-fllj]. This phenomenon can be obtained 
in pure form if the imaginary part of etotai is absent. 
If it were so, then the condition Re£ to tai < would 
lead to the unit value of the reflection coefficient, i.e. 
KV^tZ- l)/( v / £total+ 1)| 2 = 1- However, as Fig. Mb) 
shows, Ime t otai increases for greater lm£. As a result, the 
reflection for large Ime to tai is far from unity and it only 
slightly differs for both branches of the bistability (see 
Fig. The difference between the branches tends to 
decrease for Reetotai and Ime to tai as well. 



IV. NONLOCAL BISTABILITY 
(PROPAGATION EFFECTS) 

Now let us consider the layer of the two-level medium 
with finite thickness. In this case we have to take into 
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FIG. 8: Dependence of reflection on the light amplitude. The 
enhancement factor (a) I = 2 + O.li, (b) I = 2 + 0.2i, (c) 
1 = 2 + 0.3i, (d) £ = 2 + 0.4i. The other parameters are the 
same as in Fig. [7] 



account the light propagation effects resulting in the so- 
called nonlocal optical bistability [ll|. To analyze prop- 
agation of monochromatic light wave in the nonlinear 
layer, we use the iteration matrix method which was de- 
scribed in detail in Ref. [|| . This method is the modifica- 
tion of the widely known transfer matrix approach. Due 
to this method, we divide the layer into sublayers with 
thicknesses much less than the light wavelength so that 
the properties of each sublayer can be treated as constant. 
Then, according to initial conditions (refractive index of 
unexcited medium) , we calculate the distribution of light 
inside the layer and, hence, obtain new values of refrac- 
tive index of each sublayer. These calculations can be 
continued many times until, at a certain iteration, the 
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distribution of light and refractive index take on the sta- 
tionary form (with a certain accuracy). 

The results of calculations of reflection and transmis- 
sion coefficients are demonstrated in Fig. [JJ and HJ In 
Fig. [JJ one can see one more confirmation of the role of 
the local-field enhancement factor I in bistability appear- 
ance. There is also multiplicity of hysteresis loops that 
can be treated as the main feature of nonlocal bistability 
Q . The number of loops depends on the layer thickness 
and the value of i [compare figures [7Jb) and (c)]. Trans- 
mission has only one loop corresponding to the most in- 
tensive reflection loop [Fig. UJd)]. 

Even small absorption of host dielectric significantly 
changes this bistable response as shown in Fig. |U As 
the imaginary part of t grows, the loops tend to stretch 
along the light amplitude axis. The number of loops and 
their width are decreasing as well. There is another im- 
portant consequence of background absorption which is 
connected with the stability of the effect. As one can 
see in Fig. [JJb) and (c), reflection at intensities between 
the loops demonstrate many spikes distributed chaoti- 
cally. These spikes correspond to the regions of miscon- 
vergence of the iteration matrix method so that we have 
to stop it at a certain iteration and take next intensity 
value. As opposed to this behavior, the regions of bista- 



bility have rather fast convergence. Perhaps, the spikes 
are due to instability by auto-oscillations mechanism re- 
ported in [T3, EH ■ Figure [5] shows that even low level of 
absorption leads to elimination of spikes, i.e. response 
of the system becomes stable. This means that the two- 
level medium with the background dielectric of quite low 
absorption can be used to obtain only slightly suppressed 
multiple bistability without auto-oscillations. 



V. CONCLUSION 

In this article we have considered semiclassically the 
general situation of the dense collection of two-level 
atoms placed into absorbing background dielectric. Of 
course, the coefficient of dielectric absorption cannot take 
on arbitrary values. We also do not take into account the 
problem of medium heating due to radiation absorption. 
Nevertheless, the analysis carried out above can be useful 
for proper selection of the media in possible experiments. 
In particular, the appropriate choice of slightly absorb- 
ing (and, hence, slightly heated) background medium can 
help to obtain the hysteresis response without such insta- 
bilities as auto-oscillations. 
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